PCR biases distort bacterial and archaeal community structure in pyrosequencing datasets

As 16S rRNA gene targeted massively parallel sequencing has become a common tool for microbial diversity investigations, numerous advances have been made to minimize the influence of sequencing and chimeric PCR artifacts through rigorous quality control measures. However, there has been little effort towards understanding the effect of multi-template PCR biases on microbial community structure. In this study, we used three bacterial and three archaeal mock communities consisting of, respectively, 33 bacterial and 24 archaeal 16S rRNA gene sequences combined in different proportions to compare the influences of (1) sequencing depth, (2) sequencing artifacts (sequencing errors and chimeric PCR artifacts), and (3) biases in multi-template PCR, towards the interpretation of community structure in pyrosequencing datasets. We also assessed the influence of each of these three variables on α- and β-diversity metrics that rely on the number of OTUs alone (richness) and those that include both membership and the relative abundance of detected OTUs (diversity). As part of this study, we redesigned bacterial and archaeal primer sets that target the V3–V5 region of the 16S rRNA gene, along with multiplexing barcodes, to permit simultaneous sequencing of PCR products from the two domains. We conclude that the benefits of deeper sequencing efforts extend beyond greater OTU detection and result in higher precision in β-diversity analyses by reducing the variability between replicate libraries, despite the presence of more sequencing artifacts. Additionally, spurious OTUs resulting from sequencing errors have a significant impact on richness or shared-richness based α- and β-diversity metrics, whereas metrics that utilize community structure (including both richness and relative abundance of OTUs) are minimally affected by spurious OTUs. However, the greatest obstacle towards accurately evaluating community structure are the errors in estimated mean relative abundance of each detected OTU due to biases associated with multi-template PCR reactions

The active microbial community more accurately reflects the anaerobic digestion process: 16S rRNA (gene) sequencing as a predictive tool

Background: Amplicon sequencing methods targeting the 16S rRNA gene have been used extensively to investigate microbial community composition and dynamics in anaerobic digestion. These methods successfully characterize amplicons but do not distinguish micro-organisms that are actually responsible for the process. In this research, the archaeal and bacterial community of 48 full-scale anaerobic digestion plants were evaluated on DNA (total community) and RNA (active community) level via 16S rRNA (gene) amplicon sequencing. Results: A significantly higher diversity on DNA compared with the RNA level was observed for archaea, but not for bacteria. Beta diversity analysis showed a significant difference in community composition between the DNA and RNA of both bacteria and archaea. This related with 25.5 and 42.3% of total OTUs for bacteria and archaea, respectively, that showed a significant difference in their DNA and RNA profiles. Similar operational parameters affected the bacterial and archaeal community, yet the differentiating effect between DNA and RNA was much stronger for archaea. Co-occurrence networks and functional prediction profiling confirmed the clear differentiation between DNA and RNA profiles. Conclusions: In conclusion, a clear difference in active (RNA) and total (DNA) community profiles was observed, implying the need for a combined approach to estimate community stability in anaerobic digestion

Probabilistic models to describe the dynamics of migrating microbial communities

In all but the most sterile environments bacteria will reside in fluid being transported through conduits and some of these will attach and grow as biofilms on the conduit walls. The concentration and diversity of bacteria in the fluid at the point of delivery will be a mix of those when it entered the conduit and those that have become entrained into the flow due to seeding from biofilms. Examples include fluids through conduits such as drinking water pipe networks, endotracheal tubes, catheters and ventilation systems. Here we present two probabilistic models to describe changes in the composition of bulk fluid microbial communities as they are transported through a conduit whilst exposed to biofilm communities. The first (discrete) model simulates absolute numbers of individual cells, whereas the other (continuous) model simulates the relative abundance of taxa in the bulk fluid. The discrete model is founded on a birth-death process whereby the community changes one individual at a time and the numbers of cells in the system can vary. The continuous model is a stochastic differential equation derived from the discrete model and can also accommodate changes in the carrying capacity of the bulk fluid. These models provide a novel Lagrangian framework to investigate and predict the dynamics of migrating microbial communities. In this paper we compare the two models, discuss their merits, possible applications and present simulation results in the context of drinking water distribution systems. Our results provide novel insight into the effects of stochastic dynamics on the composition of non-stationary microbial communities that are exposed to biofilms and provides a new avenue for modelling microbial dynamics in systems where fluids are being transported

The active microbial community more accurately reflects the anaerobic digestion process : 16S rRNA (gene) sequencing as a predictive tool

Amplicon sequencing methods targeting the 16S rRNA gene have been used extensively to investigate microbial community composition and dynamics in anaerobic digestion. These methods successfully characterise amplicons, but do not distinguish micro-organisms that are actually responsible for the process. In this research, the archaeal and bacterial community of 48 full-scale anaerobic digestion plants was evaluated on DNA (total community) and RNA (active community) level via 16S rRNA (gene) amplicon sequencing. A significantly higher diversity on DNA compared with the RNA level was observed for archaea, but not for bacteria. Beta diversity analysis showed a significant difference in community composition between the DNA and RNA of both bacteria and archaea. This related with 25.5 and 42.3% of total OTUs for bacteria and archaea, respectively, that showed a significant difference in their DNA and RNA profiles. Similar operational parameters affected the bacterial and archaeal community, yet, the differentiating effect between DNA and RNA was much stronger for archaea. In conclusion, a clear difference in active (RNA) and total (DNA) community profiles was observed, implying the need for a combined approach to estimate community stability in anaerobic digestion

Anaerobic microbial community response to methanogenic inhibitors 2‐bromoethanesulfonate and propynoic acid

Methanogenic inhibitors are often used to study methanogenesis in complex microbial communities or inhibit methanogens in the gastrointestinal tract of livestock. However, the resulting structural and functional changes in archaeal and bacterial communities are poorly understood. We characterized microbial community structure and activity in mesocosms seeded with cow dung and municipal wastewater treatment plant anaerobic digester sludge after exposure to two methanogenic inhibitors, 2‐bromoethanesulfonate (BES) and propynoic acid (PA). Methane production was reduced by 89% (0.5 mmol/L BES), 100% (10 mmol/LBES), 24% (0.1 mmol/LPA), and 95% (10 mmol/LPA). Using modified primers targeting the methyl‐coenzyme M reductase (mcrA) gene, changes in mcrA gene expression were found to correspond with changes in methane production and the relative activity of methanogens. Methanogenic activity was determined by the relative abundance of methanogen 16S rRNA cDNA as a percentage of the total community 16S rRNA cDNA. Overall, methanogenic activity was lower when mesocosms were exposed to higher concentrations of both inhibitors, and aceticlastic methanogens were inhibited to a greater extent than hydrogenotrophic methanogens. Syntrophic bacterial activity, measured by 16S rRNA cDNA, was also reduced following exposure to both inhibitors, but the overall structure of the active bacterial community was not significantly affected.This manuscript reports a comprehensive approach to characterizing the effects of commonly used methanogenesis inhibitors on an anaerobic microbial community. We use mock and environmental communities and target two genes using DNA‐ and RNA‐based methods. Results from Illumina sequencing of the 16S rRNA gene, 16S rRNA cDNA, mcrA gene, and mcrA transcript cDNA highlight shifts in both methanogenic archaeal activity and syntrophic bacterial activity.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/134127/1/mbo3349.pdfhttp://deepblue.lib.umich.edu/bitstream/2027.42/134127/2/mbo3349_am.pd

Diverse manganese(II)‐oxidizing bacteria are prevalent in drinking water systems

Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/136363/1/emi412508_am.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/136363/2/emi412508.pd

Metagenomic evidence for the Presence of Comammox Nitrospira-Like Bacteria in a Drinking Water System

We report metagenomic evidence for the presence of a Nitrospira-like organism with the metabolic potential to perform the complete oxidation of ammonia to nitrate (i.e., it is a complete ammonia oxidizer [comammox]) in a drinking water system. This metagenome bin was discovered through shotgun DNA sequencing of samples from biologically active filters at the drinking water treatment plant in Ann Arbor, MI. Ribosomal proteins, 16S rRNA, and nxrA gene analyses confirmed that this genome is related to Nitrospira-like nitrite-oxidizing bacteria. The presence of the full suite of ammonia oxidation genes, including ammonia monooxygenase and hydroxylamine dehydrogenase, on a single ungapped scaffold within this metagenome bin suggests the presence of recently discovered comammox potential. Evaluations based on coverage and k-mer frequency distribution, use of two different genome-binning approaches, and nucleic acid and protein similarity analyses support the presence of this scaffold within the Nitrospira metagenome bin. The amoA gene found in this metagenome bin is divergent from those of canonical ammonia and methane oxidizers and clusters closely with the unusual amoA gene of comammox Nitrospira. This finding suggests that previously reported imbalances in abundances of nitrite- and ammonia-oxidizing bacteria/archaea may likely be explained by the capacity of Nitrospira-like organisms to completely oxidize ammonia. This finding might have significant implications for our understanding of microbially mediated nitrogen transformations in engineered and natural systems

Gradual Recovery of Building Plumbing-Associated Microbial Communities after Extended Periods of Altered Water Demand during the COVID-19 Pandemic

COVID-19 pandemic-related building restrictions heightened drinking water microbiological safety concerns post-reopening due to the unprecedented nature of commercial building closures. Starting with phased reopening (i.e., June 2020), we sampled drinking water for 6 months from three commercial buildings with reduced water usage and four occupied residential households. Samples were analyzed using flow cytometry and full-length 16S rRNA gene sequencing along with comprehensive water chemistry characterization. Prolonged building closures resulted in 10-fold higher microbial cell counts in the commercial buildings [(2.95 ± 3.67) × 105 cells mL-1] than in residential households [(1.11 ± 0.58) × 104 cells mL-1] with majority intact cells. While flushing reduced cell counts and increased disinfection residuals, microbial communities in commercial buildings remained distinct from those in residential households on the basis of flow cytometric fingerprinting [Bray-Curtis dissimilarity (dBC) = 0.33 ± 0.07] and 16S rRNA gene sequencing (dBC = 0.72 ± 0.20). An increase in water demand post-reopening resulted in gradual convergence in microbial communities in water samples collected from commercial buildings and residential households. Overall, we find that the gradual recovery of water demand played a key role in the recovery of building plumbing-associated microbial communities as compared to short-term flushing after extended periods of reduced water demand

Drift dynamics in microbial communities and the effective community size

The structure and diversity of all open microbial communities are shaped by individual births, deaths, speciation and immigration events; the precise timings of these events are unknowable and unpredictable. This randomness is manifest as ecological drift in the population dynamics, the importance of which has been a source of debate for decades. There are theoretical reasons to suppose that drift would be imperceptible in large microbial communities, but this is at odds with circumstantial evidence that effects can be seen even in huge, complex communities. To resolve this dichotomy we need to observe dynamics in simple systems where key parameters, like migration, birth and death rates can be directly measured. We monitored the dynamics in the abundance of two genetically modified strains of Escherichia coli, with tuneable growth characteristics, that were mixed and continually fed into 10 identical chemostats. We demonstrated that the effects of demographic (non‐environmental) stochasticity are very apparent in the dynamics. However, they do not conform to the most parsimonious and commonly applied mathematical models, where each stochastic event is independent. For these simple models to reproduce the observed dynamics we need to invoke an “effective community size”, which is smaller than the census community size